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We present a study of the effects of improvement on light-quark physics using the Fermilab formalism. The 
calculations were performed at three different lattice spacings, with the SW action and traditional Wilson action 
(both tadpole improved). We find that O(a) effects for the decay constant f w and quark mass m q can be 
successfully removed using the tadpole-improved SW action in the light-quark regime. 



1. Introduction 

In this project, we have made a systematic 
comparison of improved and unimproved actions 
by studying the lattice-spacing dependence of 
the simplest quantities in light-quark physics: 
the pion decay constant, f w , and the (isospin- 
averaged) light-quark mass, m q . 

The first step in improving the Wilson action 
for QCD is to remove the leading 0(a) errors 
from Wilson fcrmions. Following Symanzik, this 
leads to the Sheikholeslami-Wohlert (SW) 0(a)- 
improved fermion action [EJ. With tadpole im- 
provement ||, this gives a QCD action which is 
essentially correct to 0(a 2 ). 

A reliable improved action for light quarks 
would allow unquenched calculations on small 
lattices, with important implications for phe- 
nomenology. 

2. Lattice Details 

We study scaling effects using three lattices, at 
f3 = 5.7, 5.9, and 6.1, which cover a range of more 
than a factor of two in a. 

On each lattice we compute Wilson and SW 
propagators, as described below. In this study we 
aim to test the effects of Symanzik improvement, 
i.e. the presence or absence of the SW clover O(a)- 
improvement term in the action. We do not in- 
vestigate the effects of tadpole improvement com- 
pared with no tadpole improvement — tadpole 
improvement is used throughout, and the same 
tadpole-improvement procedure is used for both 
actions. 

The full lattice details are shown in Table [j]. 
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Lattice details (n/ = 0) 




0=6.1 


/3=5.9 


0=5.7 


L 3 x T 


24 3 x 48 


16 3 x 32 


12 3 x 24 


a- 1 


2.43 GeV 


1.78 GeV 


1.15 GeV 


JJihys 


2.0 fm 


1.8 fm 


2.1 fm 


C 


1.46* 


1.50* 


1.57 


configs 








Wilson 


100 


100 


200 


SW 


100 


150 


200 


# light k's 


5 


5 


5 



*For historical reasons, the (3 = 6.1 and 5.9 results were 
obtained with the mean-field value c = 1.4. The differ- 
ences are negligible and do not modify the conclusions. 



3. Improvement Procedure 

We use a tadpole-improved SW fermion action, 

S = S Wllson + (i/2) C^G-F^, (1) 

where the coefficient of the clover term, c, is set 
non-perturbatively, c = l/(u ) 3 . The average 
link, uq, is determined from the plaquette, 

Mo = ((l/3)Trcy 1/4 (2) 

to give the values of c shown in Table |l|. The 
original tree-level SW action used c = 1. 

For the fermions, we take the tadpole-improved 
field ip, 

%j} = y/l- 3K/4K c „ t -0 (3) 
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which has the correct normalization in the chiral 
limit, tp = "0/2. 

In calculating matrix elements with the SW ac- 
tion an improved operator must be used ||. In 
the formalism of Heatlie et al this is equivalent 
to rotating the quark fields, 

[l-( z p-{l-z)mo)/2] i\>. (4) 

Choosing z — 0, the rotation becomes, 

ip -» [1 + m /2] </> (5) 

and its effects will vanish in the chiral limit, where 

?7lo = 0. 

With the Fermilab formalism the correspond- 
ing rotation is, 



1p -» ^/T+ 



m 



1 + d x D ■ 7 



(6) 



but its effects also vanish in the chiral limit (since 
d\ ~ mo for small too). For convenience we 
used this trick to obtain our decay-constant re- 
sults from the ordinary unrotated fields. 

In chiral extrapolations we use linear fits in the 
lattice pole mass, 



m q = In (1 + (1/2k - l/2K crit )) 



(7) 



to model the K-dependence. This definition au- 
tomatically resums the tree-level large-m correc- 
tions in the quark mass. 

For the Wilson action we use exactly the same 
tadpole-improvement procedure for the fermions, 
but set c = 0. 

4. Setting the Scale 

In a systematic study of 0(a) effects it is im- 
portant to choose a scale which is as free from 
a-dependence as possible, so that other quanti- 
ties are not affected by a common variation in 
the scale. 

We set the scale using the spin-averaged IP — 
IS splitting from charmonium jj]. The spin- 
averaged splitting is found to be insensitive to 
0(a) corrections from the clover term, and the 
resulting scale is free from 0(a) effects. Using 
this scale, the a-dependence of other quantities 
can be determined directly, up to the remaining 
0(a 2 ) systematic errors in the scale. 



The advantage of using charmonium to set the 
scale is that it is possible to separate the a- 
dependence of the quantities under study from 
the a-dependence of the scale (known to be good 
to 0(a 2 )). Other possible scales, such as m p , have 
the disadvantage that the scale itself can contain 
an 0(a) effect. 

5. Multistate Fitting 

When studying small 0(a) effects it is im- 
portant to keep systematic errors under control. 
These errors may have an a-dependence which 
conceals the true scaling behavior of the quan- 
tities under study. In particular, excited-state 
contamination can introduce large 0(a) effects, 
purely from the systematic variations in a poor 
fitting procedure. 

We use a sophisticated multistate fitting pro- 
gram to eliminate excited-state contamination. A 
set of optimised smearing functions is available 
on each lattice, created from a study of coulomb- 
gauge pion wavefunctions, typically for IS and 
2S sources. The fitting is performed using a ma- 
trix of correlators, 6-6 (point source), 1S-1S, 25- 
2S, . . . , and a choice of one-exponential or 

two-exponential fits. We check for consistency be- 
tween the different fitting methods, to ensure that 
excited-state contamination has been removed. 

6. Decay Constant 

Previous studies of f n using Wilson fermions 
have all shown a negligible dependence on lattice- 
spacing. In order to determine whether this is 
due to an intrinsically small 0(a) correction for 
/ w from the action, or a cancellation by 0(a) ef- 
fects from excited-state contamination, our anal- 
ysis was carried out using multistate fits. 

Fermilab multistate results for f n using the 
tadpole-improved Wilson and tadpole-improved 
SW actions are shown in Fig. ([!]), with the phys- 
ical value /,r = 132 MeV. The lattice values are 
also extrapolated to the continuum limit using a 
simple linear fit in a. 

The Wilson results are renormalized using the 
tadpole-improved Za of Lepage and Macken- 
zie H . The tadpole- improved value of Za for the 
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Figure 1 . Comparison of S W and Wilson / w . The 
lines show the effect of varying q* from 1/a (lower 
line) to it /a (upper line). 



Figure 2. Comparison of SW and Wilson m q . The 
lines show the effect of varying q* from 1/a (lower 
line at large a) to it /a (upper line at large a). 



local current with the SW action was obtained 
from the work of Borrelli et al Q , by subtracting 
the tadpole contribution from their intermediate 
results, 

Z A (q*) = 1 - 0.31a v (q*) (Wilson) (8) 
= 1 - 0.61a v (q*) (SW). (9) 

The exact values of q* remain to be determined by 
analytic calculations. For the moment, we allow 
q* to vary over a reasonable range, from 1/a to 
it /a, and regard this as a systematic error. 

The additional systematic errors from the IP — 
IS charmonium scale are expected to be 0(p 2 a 2 ), 
which at /3 = 5.7 might be 5-10%. 

Allowing for these systematic errors, the mul- 
tistate fits confirm the behavior of the existing 
Wilson data, and show a small lattice-spacing de- 
pendence. This indicates that the coefficient of 
the 0(a) correction to Wilson is intrinsically 
small. 

The SW results are consistent with a small or 
no 0(a) effect within systematic errors. The ad- 
dition of the O(a) correction term to the Wil- 



son action does not increase a-dependence of the 
results when the a-dependence is already small. 
This is a requirement of improvement. 

Even allowing for the uncertainty in q* , the lin- 
ear extrapolations of Wilson and SW f v to the 
continuum limit do not meet. This is a possible 
indication of the perturbative (g 4 ) or quadratic 
(a 2 ) corrections to the naive linear scaling law. 

7. Quark Mass 

We extract the light-quark mass from the pion 
mass chiral extrapolation (by determining the 
quark mass that gives the physical pion mass), 

ml = Am q . (10) 

In Gupta's Lattice '94 review of quark masses 
a compilation of existing Wilson data showed a 
large variation of m q with lattice-spacing, and a 
significant discrepancy between Wilson and Stag- 
gered results H . 

Again, in order to determine whether this is 
due to systematic 0(a) effects from excited-state 
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contamination, our analysis was carried out using 
multistate fits. 

The Fermilab multistate results for m q using 
tadpole-improved Wilson and tadpole-improved 
SW actions are shown in Fig.(||). The lattice m q 
values are converted to a common scale of (3 = 5.9 
using the 1-loop relation, 

m(a 5 . 9 ) = (1 - (2/w)a(q*) ln(a/a 5 . 9 )) m(a). (11) 

This removes the lattice-spacing dependence due 
to the running of the mass. We keep the masses 
in the lattice scheme. An additional constant fac- 
tor, Zj-j-g, would be needed to obtain the conven- 
tional M S value. Since we are only interested in 
the scaling properties this overall constant is not 
important here. 

The masses are again extrapolated to the con- 
tinuum limit using a simple linear fit in a. In 
this case the Wilson and SW extrapolations are 
in good agreement. 

Multistate fitting confirms the presence of a 
large 0(a) effect in Wilson m q . The SW re- 
sults show a small O(a) effect, compatible with 
zero within systematic errors. This indicates that 
the discrepancy between the existing Wilson and 
Staggered results is probably due to a large 0(a) 
error from the Wilson action. The Staggered ac- 
tion, which is also correct to 0(a 2 ), should be in 
agreement with the tadpole-improved SW action. 

8. Conclusions 

For the decay constant, multistate fitting con- 
firms the small O(a) effects found in Wilson 
The O(a) effects for the improved SW f„ are also 
small, within systematic errors from q* and p 2 a 2 
corrections to the charmonium scale. There is a 
difference between the continuum extrapolations 
of SW and Wilson / w which suggests perturbative 
or quadratic corrections to naive linear scaling. 

For the quark mass, multistate fitting confirms 
the large O(a) effects found in Wilson m q . The 
0(a) effects for the SW m q are consistent with 
zero, within systematic errors. The improved ac- 
tion removes the a-dependence of m q even down 
to (3 = 5.7, when using the clover term with tad- 
pole improvement. 

Further simulations at lower [3 values are desir- 



able, to determine the scale at which 0(a 2 ) errors 
become important. At this stage there is no clear 
evidence of 0(a 2 ) effects, even at [3 — 5.7. How- 
ever, for larger lattice spacings it is essential to 
first obtain the correct q* values, since the sys- 
tematic error due to q* increases as the perturba- 
tive corrections increase. 

The costs involved in using the SW action 
are small compared with the benefits of reducing 
0(a) errors. On the canopy/ ACPMAPS platform 
the addition of the clover term typically increases 
the runtime of the quark inversion by approxi- 
mately 30%. 

In summary, when Wilson results already show 
small 0(a) effects, improvement does not make 
things worse, and when Wilson results show large 
0(a) effects, improvement reduces them. Thus, 
use of the tadpole-improved SW action is recom- 
mended. 
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